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We consider a family of random walks killed at the boundary of the Weyl chamber 
of the dual of Sp(4), which in addition satisfies the following property: for any n ^ 3, 
there is in this family a walk associated with a reflection group of order 2n. Moreover, 
the case n = 4 corresponds to a process which appears naturally by studying quantum 
random walks on the dual of Sp(4). For all the processes belonging to this family, we 
find the exact asymptotic of the Green functions along all infinite paths of states as well 
Qu( ■ as that of the absorption probabilities along the boundaries. 

Resume 

Dans cet article, nous considerons une famille de marches aleatoires tuees au bord 
de la chambre de Weyl du dual de Sp(4), qui verifie en outre la propriete suivante : pour 
tout n ^ 3, il y a, dans cette famille, une marche ayant un groupe de reflexions d'ordre 
2n. De plus, le cas n — 4 correspond a un processus bien connu apparaissant lors de 
1'etude des marches aleatoires quantiques sur le dual de Sp(4). Pour tous les processus de 
cette famille, nous trouvons l'asymptotique exacte des fonctions de Green selon toutes 
les trajectoires, ainsi que l'asymptotique des probabilites d'absorption sur le bord. 
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1 Introduction and main results 

Appearing in several distinct domains, random walks conditioned on staying in cones of Z 
attract more and more attention from the mathematical community. Historically, important 
examples are the so-called non-colliding random walks. These are the processes (Zi, . . . , Zj) 
composed of d independent and identically distributed random walks conditioned on never 
leaving the Weyl chamber {z £ M d : Z\ < ■ ■ ■ < Zd}- They first appeared in the eigenvalues 
description of important matrix- valued stochastic processes, see [11], and are recently again 
very much studied, see |12} [8j [21] and the references therein. Another important area where 
processes conditioned on never leaving cones of Z ' appear is that of quantum random walks, 
see e.g. PH]. 

A usual way to condition random processes on staying in cones consists in using Doob 
/i-transforms. These are functions which are harmonic, positive inside of the cone and equal 
to zero on its boundary — or equivalently harmonic and positive for the underlying killed 
processes. It is therefore natural to be interested in finding all positive harmonic functions 
for processes in cones of Z d killed at the boundary, and more generally to compute the Martin 
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compactification of such processes, that can e.g. be obtained from the exact asymptotic of 
the Green functions. 

We briefly recall [10] that for a transient Markov chain with state space E, the Martin 
compactification of E is the smallest compactification E of E for which the Martin kernels 
y h4 ky = Gy/ Gy° extend continuously — by Gy we mean the Green functions and we denote 
by xq a reference state. E\E is usually called the full Martin boundary. For a £ E, x h4 
is clearly superharmonic; then d m E = {a € E \ E : x \— > is minimal harmonic} is called 
the minimal Martin boundary — a harmonic function h is said minimal if ^ h ^ h with h 
harmonic implies h = ch for some constant c. Then, every superharmonic (resp. harmonic) 
function h can be written as h(x) = Jgkyfi(dy) (resp. h(x) = f d E kyfi(dy)), where \i is 
some finite measure, uniquely characterized in the second case above. 

In this context, the case of walks in cones of Z d spatially homogeneous in the interior, 
with non-zero drift and killed at the boundary has held a great and fruitful deal of attention. 

For random walks on weight lattices in Weyl chambers of Lie groups, Biane [4j finds the 
minimal Martin boundary thanks to Choquet-Deny theory. Collins [7] obtains their Martin 
compactification. 

In [23] we give a more detailed analysis for a certain class of walks in dimension d = 2. 
These are the random walks killed at the boundary of Z+, with non-zero jump probabilities 
to the eight nearest neighbors and having in addition a positive mean drift. The asymptotic 
of the Green functions along all infinite paths of states as well as that of the probabilities of 
absorption along the axes are computed for the walks in this class. However, the methods of 
complex analysis used in [23] apply in dimension d = 2 only. 

Ignatiouk- Robert [15, 14J, then Ignatiouk-Robert and Loree [16J find the Martin compac- 
tification of the random walks in Z + x Z^ -1 and lA (d #s 2), with non-zero drift and killed at 
the boundary. They make very general assumptions on the jump probabilities. The approach 
used there, based on large deviation techniques and Harnack inequalities, seems not to be 
powerful for studying the asymptotic of the Green functions. Furthermore, having a non- 
zero drift is an essential hypothesis in [151 HI El- Last but not least, the results of [H] in 
the case d ^ 3 are conditioned by the fact of been able: "to identify the positive harmonic 
functions of a random walk on Z^ which has zero mean and is killed at the first exit from ; 
unfortunately, for d ^ 2, there are no general results in this domain" (see page 5 of [14]). 

As may this open problem suggest, the results and methods dealing with the asymptotic 
of Green functions or even with the Martin compactifiction for random walks in domains of 
Z d with drift zero and killed at the boundary are actually scarce, even for d = 2. 

In [25] (resp. [22]), approximations of the Green functions for simple random walks killed 
at the boundary of balls of 7h d (resp. of certain more general sets of Z 2 ) are computed by 
comparison with Brownian motion. Namely, the Green functions G{B)y of Brownian motion 
(resp. random walk) killed at the boundary B are related to the potential kernels (or the 
Green functions if d ^ 3) a* of the non-killed Brownian motion (resp. random walk), e.g. 
via the "balayage formula" , see Chapter 4 of [25J . These identities have the same form for 
Brownian motion and random walk [25j: 

G(Bf y = -a* + E x [a S y TB ] + F(BT, 

where S is the process, tb is the hitting time of the boundary B and where the rest F{B) X 
can be expressed in terms of tb and x, indeed see Chapter 4 of [25]. They are then compared 
term by term. For the comparison of potential kernels a x of the non-killed processes, classical 
formulas may be used, like that stated in Chapter 4 of [25J; for the comparison of positions of 
these processes at time of absorption, strong approximations of random walks by Brownian 
motion [HI [20] as well as Beurling estimates [T8J [24] are usually exploited. 

This approach in particular requires rather precise estimates of the hitting time tb- It 
seems therefore difficult to use for models of walks in cones where tb is painful to analyze. 
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This is for example the case of the random walks in the half-plane considered in [29\ [30] and 
of the random walks in the quarter-plane we shall study in this paper (for which Chapter F 
of [27] illustrates the complexity of hitting times). Remark[8]in Subsection ll,3l specifies other 
reasons why this method via comparison with Brownian motion seems not to lead to enough 
satisfactory results in the analysis of the random walks we shall consider here. 

For random walks in the half-plane, Uchiyama [29, 30j finds the asymptotic of the Green 
functions from their trigonometric representations. 

If the domain is a quarter-plane, the simplest case is the cartesian product of two killed 
one-dimensional simple random walks with mean zero. The Martin boundary then happens 
to be trivial [26]. Moreover, the exact asymptotic of the Green functions for these processes 
is computed in Chapter D of [27] , 

From a Lie group theory point of view, the previous case corresponds to the group 
product SU(2)xSU(2), which is associated with a reducible rank-2 root system, see [5]. We 
are then interested in the classical random walks that can be obtained from the construction 
made by Biane in [3] — namely, by restriction of the quantum walks, initially defined on 
non-commutative von Neumann algebras, to commutative subalgebras — starting from Lie 
groups associated with irreducible rank-2 roots systems, namely SU(3) and Sp(4). These are 
random walks on the lattices of Figure [IJ spatially homogeneous in the interior and absorbed 
at the boundary. The one on the left, associated with SU(3), has three jump probabilities 
equal to 1/3; that on the right, related to Sp(4), has the same four jump probabilities 1/4. 




Figure 1: Random walks in the Weyl chamber of the duals of SU(3) and Sp(4) 

By obvious transformations of these lattices it is immediate that both are killed nearest 
neighbors random walks in the quarter-plane as below. 




Figure 2: Random walks of Figure [T] can be viewed as random walks in 

For the random walk associated with SU(3) (and besides also for its multi-dimensional 
analogues SU((i)), Biane [3] computes the asymptotic of the Green functions along all paths 
of states except for the ones approaching the axes. To complete the latter results for these 
particular paths, the complex analysis methods of [28] recently turned out to be fruitful. 
Furthermore, in [28] we compute the exact asymptotic of the Green functions along all paths 
and also the asymptotic of the absorption probabilities for a whole class of nearest neighbors 
random walks killed at the boundary of Z^_. This class, including the walk associated with 
SU(3), is characterized by the fact that each of its elements has a harmonic function of 
order three, i.e. (io, jo) i-» iojo(io + cr/o + P) for some arbitrary a and f3. Similar results for 
the simpler class of walks admitting the harmonic function of order two (io> Jo) l— Wo are 
derived in [27] . 

The methods of [28, 27J heavily rely on the analytic approach for the random walks in the 
quarter-plane developed by Fayolle, Iasnogorodski and Malyshev [13]. An important notion 
involved there is that of the group of the random walk, which is a group of automorphisms of 
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an algebraic curve; it is here properly defined in Subsection ll.il In this context our previous 
works [28\ [27] treat the classes of walks with groups of the smallest orders four and six. 

The main subject of this paper is to obtain, by this approach, the exact asymptotic of 
the Green functions along all paths as well as that of the absorption probabilities along the 
axes for a certain class of random walks in absorbed at the axes, with zero drift and 
containing the walk associated with Sp(4) — drawn on the right of Figured This class will 
include walks with groups of all finite orders 2 re for n ^ 3. We now define this class and 
specify the notion of group of the random walk. 

1.1 Random walks under consideration and group of the walk 

Consider the random walk (X(k),Y(k))k^o spatially homogeneous inside of the quarter- 
plane Z 2 and such that if p id = F[(X(k + l),Y(k + 1)) = (i + i,j + j) | (X(k),Y(k)) = 
(io, jo)], then: 

(Hi) pi )0 + pi _i + p-i,o + p-i,i = l, pi,o = p-1,0, Pi ,-i = 

(H2) {(i,0) : i ^ 1} U {(0, j) : j > 1} is absorbing. 

Let us also define the polynomial Q (which is just a simple transformation of the transition 
probabilities generating function) by: 

Q(x,y) = xy[pi i0 x + p-ifl/x + pi-ix/y + p-i t iy/x - l] . (1) 

If Q(x,y) = 0, then with ([lj it is immediate that Q(£(x,y)) = and Q(rj(x,y)) = 0, where 

y) = ( x, — V ri(x,y) = f p-^y+p-^ y^y 

V y J V Pi,oy + pi -i x 

The group of the walk is then W = (^,t(), the group of automorphisms of the algebraic 
curve {(x,y) G (Cu{oo}) 2 : Q{x,y) = 0} generated by £ and rj. Its order is always even and 
larger than or equal to four. It is already known [6] that if pi,o = P-i,o = 1/4 and pi-i = 
= 1/4 then W has order eight. 

More generally, we prove in Remark[TT]that the group W is finite if and only if there exists 
some rational number r such that p\$ = p-ip = sin(r7r) 2 /2 and pi-% = = cos(r7r) 2 /2. 

As illustrated by the works already mentioned \1'6\ El I2"8] [2T] the notion of finite group 
is nowadays extensively studied, notably because it often leads to worthwhile results. Let 
&2n be the class of random walks satisfying (HI), (H2) and with a group of order 2n. We 
show in Subsection 12.21 that the random walk under hypotheses (HI), (H2) and (H3), where 

(H3) pi j0 = P-i,o = sin(7r/n) 2 /2 and p 1 _ x = p_ 1)X = cos(7r/n) 2 /2, 

belongs to &2n- 

In this article we study the class being made up of the union for n ^ 3 of the random walks 
satisfying (HI), (H2) and (H3). This class contains one — and only one — representative of 
^2n for any n ^ 3. Precisely, for all walks in this class we compute the exact asymptotic of 
the Green functions along all paths and that of the absorption probabilities along the axes. 
This is the first result of that kind for random walks with zero drift and groups of all finite 
orders, up to our knowledge. 

In the particular case n = 3, the process coincides with that represented on the left of 
Figure 2 in [28J for a = 2. In the case n = 4, this is the walk in the Weyl chamber of Sp(4) 
with jump probabilities 1/4 studied by Biane [3], see Figures H] and [2] 

The hypothesis that n is integer is technical. Indeed, independently of this assumption, 
the approach [13] — we shall use here — always yields explicit expressions for the Green 
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functions, notably in terms of solutions to boundary value problems of Riemann-Hilbert 
type. In the general case, these formulations are so complex that we are not able to obtain 
their asymptotic. However, they may admit a nice simplification as a closed expression (and 
then in terms of the orbit under the group of a simple function); this actually happens if 
the walk admits a finite group, and if in addition some technical assumption holds — related 
to fundamental domains, see Subsection 12.21 As it will be properly showed in Remark II 11 
in the case of the walks satisfying to (HI) this exactly implies (H3). 



1.2 Main results 

Here and throughout, (X, Y) = (X(k),Y(k))k^o denotes the process defined by (HI), (H2) 
and (H3). Our first result deals with the asymptotic of the Green functions, properly defined 
by 



/-~l*0 JO TI7 



Yl 1 {(x(fc),y(fe))=(i,i)} 



(2) 



Let f n be the function defined in f|23|) : in Section [3] we shall write it explicitly and we shall 
prove that it is harmonic for (X, Y), positive inside of and equal to zero on the boundary. 

Theorem 1. The Green functions (0) admit the following asymptotic as i + j — >• oo and 
j/i ->tan(7), 7 € [0,tt/2]; 



^o jo 



[n 



X)l sin (narctan [ j^jji tan(7r/n)]) 



/»(*>.*>)—> — (3) 



vr 4" sin (2vr/n) ' [ CO s(vr/n) 2 (i 2 + 2ij) + j 2 ]™ /2 



Remark 2. Let N n {j/i) = sin (narctan [ \J.ju tan(7r /rt)] ) be the quantity appearing in the 
asymptotic (3\). Let also 7 be in [0, ir/2] and suppose that j/i goes to tan (7). 

If "f S]0, 7r/2[, then N n (j/i) goes to iV n (tan(7)), which belongs to ]0, oo[. 

If 1 = or 7 = 7r /2 ; then N n (j/i) goes to 0. More precisely, N n (j/i) = ntan(7r/ra)[j/i + 
O(jfi) 2 } 2/7 = and N n (j/i) = (nsin(2vr/n)/2)[i/j + 0(i/j) 2 } ifj = vr/2. 

Theorem Q] has the following immediate [10] consequence. 

Corollary 3. The Martin compactification is the one-point compactification. 

This paper therefore gives a partial answer, for d = 2, to the open problem highlighted 
by Ignatiouk- Robert in p3], since Corollary [3] implies that up to the positive multiplicative 
constants, there is only one positive harmonic function for (X,Y). 

Theorem Q] also has a consequence on the absorption probabilities P( ioJ - )[(X, Y) is killed 
at (i,0)] and P^ t j Q \[(X, Y) is killed at (0,j)]. Indeed, by using the obvious equalities 



F [io>jo) [(X,Y) is killed at (i, 0)] = Pi.-iG^, 
F [ioJo) [(X,Y) is killed at (0, j)] = P-i,iG^ x +p_ li0 Gfj 7 ' 

as well as Theorem [1] and Remark [21 we come to the following result. 

Corollary 4. The absorption probabilities admit the following asymptotic as i,j —■ 00: 



\i ,j )i( x > Y ) is kiUed at (*.0)] ~ Z[7-^t™/» (*o,Jo) 



1 n\ I 

2tt [4cos(7r/n)] n/n ^°' Jo V 
In!., ,1 



2 n /±n Jn ^ U,JU ' „-n+l • 
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1.3 Harmonic functions and link with Brownian motion 

Let us now have a closer look at the harmonic function f n that governs the asymptotic ([3]) 
of the Green functions (P) . All the results of Subsection 11.31 are proven in Section [3l 

Proposition 5. 

(i) f n is a real polynomial in the variables io,jo °f degree exactly n; 

(ii) f n is a harmonic function for the process (X,Y); 
(hi) f n (io,0) = fn(OJo) = for all integers i and j ; 

(iv) Ifi ,jo > then f n (io,jo) > 0. 

The explicit formulation of f n for general values of n — that we shall obtain in Section [3] — 
being quite complex, here we just give the following three examples: 

/3(»o,io) = 24-3 1 / 2 -wo(io + 2j ), 
fSoih) = (256/3) •ioio(«o + 2io)(io+io) ! 

MioJo) = (288/5)3 1 / 2 • wo(^o + 2io)(io+io)((*o + 2io/3)(io+4jo/3) + 10/9). 

Note that for n = 3 and n = 4, f n is a homogeneous function of degree n (i.e. / n (Ax, Ay) = 
X n f n (x,y)), while is not. In fact the next result holds. 

Proposition 6. For any n 5, f n is not homogeneous. 

Let us conclude the introduction by outlining the link of f n with the harmonic functions 
of Brownian motion. Let 

&(x,y) = ((x + y)/sin(vr/n),y/cos(7r/n)). (4) 

Then the random walk 4>(X, Y) has an identity covariance and takes its values in the cone 
A(0,7r/n) = {£exp(z6l) :0^t^co,0^6l^ vr/n}. 





Figure 3: On the left, the walk (X, Y); on the right, the walk <f>(X, Y), with <p defined in Q 

4>(X, Y) therefore lies in the domain of attraction of the standard Brownian motion killed 
at the boundary of A(0,7r/n). 

For the Brownian motion, it is well-known that there is only one harmonic function h 
positive inside of a given cone and vanishing on the boundary: it is called the reduite [1] 
of the cone. It happens to be homogeneous pQ. When the cone is A(0,7r/ra), the reduite is 
equal to h(pexp(i9)) = p n sm(n9). 

Moreover, the asymptotic of the Green functions of the Brownian motion killed at the 
boundary of A(0, n/n) can be obtained from [2] and is equal to: 

^ P exp(i0) 2 sin(nr/) 

G rexpU) ~ ^Mpexp^))— r -> CO. 

Proposition 7. Let <p be defined in Q). Up to a multiplicative constant, the homogeneous 
function h(<f>(io, jo)) equals the dominant term of the non-homogeneous harmonic function 
of the random walk, i.e. 

fn(k,jo) = h((/)(io,jo))[l +o(l)], i ,j -> oo. 
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Proposition [7] will follow from our results by a direct computation, see Section [3l Let us 
also note, as in [29], that this proposition is in accordance with Donsker's invariant principle. 

Remark 8. Proposition [?] also entails that the comparison approach of random walks by 
Brownian motion sketched in the first part of the introduction can neither give the approxi- 
mation of the Green functions with the precision of Theorem{l\nor even to specify the unique 
harmonic function for our class of random walks (HI), (H2) and (US). 

Furthermore, from the asymptotic results (dp we notice that the Green functions can 
tend to zero arbitrarily fast whenever the order of the group is taken high enough, while the 
Beurling estimates U8[ relating the random walk to the Brownian motion typically have 
a polynomial precision — which in addition depends only on the dimension. 

The rest of the paper is organized as follows. In Section[2]we find explicitly the absorption 
probabilities and the Green functions ([2]). In Section [3j we then study precisely the harmonic 
function f n . Finally in Section [J] we prove Theorem [TJ 

2 Expression of the absorption probabilities and of the Green 
functions 

2.1 A functional equation between the generating functions 

Subsection 12. II consists in preparatory results and is inspired by the book [13]. Define 

G^(x,y) = ^Gf/s'-V" 1 , 

h io,jo {x) = £> ( . oJo) [(^,y) is killed at (*,0)K, (5) 
%ioJo {y) = J- W (loJo) [(X, Y) is killed at (0, j)]y* 

m 

the generating functions of the Green functions ([2]) and of the absorption probabilities. With 
these notations, we can state the following functional equation: 

Q (x, y) G lo ' jo (x, y) = h io ' jo (x) + h io ' jo (y) - x io y jo , (6) 

Q being defined in ([I]). A priori, Equation ([6]) has a meaning in {(x,y) € C 2 : \x\ < 1, \y\ < 
1}. The proof of ([6]) is obtained exactly as in Subsection 2.1 of [23j . 

When no ambiguity on the initial state can arise, we will drop the index io,jo and we 
will write G id ,G{x,y),h{x),h{y) for Gfj jo , G io 'i°(x, y), h io ^°(x), h l °^°{y). 

Let us now have a look to the algebraic curve {(x,y) € (C U {oo}) 2 : Q(x,y) = 0}, that 
we note J2 for the sake of briefness. Start by writing the polynomial (pQ) alternatively 

Q (x, y) = a (x) y 2 + b(x)y + c{x) = a (y) x 2 + b (y) x + c{y) , (7) 

where a(x) =pi-i, b(x) = piflX 2 - x +pi )0 , c(x) = p 1 _ix 2 and a(y) =p 1 gy+pi- 1 , b(y) = 
—y, c(y) = pi-iy 2 + Pi,oy- Set also d(x) = b(x) 2 — 4a(x)c(x) and d(y) = b(y) 2 — 4a(y)c(y). 
We have 

d(x) =pl {x - l) 2 (x 2 + 2x(l - l/p 1)0 ) + 1), d(y) = -4pi fl pi,- iy (y - if. (8) 

The polynomial d has manifestly a double root at 1 and two simple roots at positive 
points, that we denote by x\ < 1 < X4. As for d, it has a double root at 1 and a simple root 
at 0. We also note y\ = and 2/4 = 00. 
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Then with ([7]) we notice that Q(x,y) = is equivalent to [b(x) + 2a(x)y] 2 = d{x) or 
to [b(y) + 2a(y)x] 2 = d(y). It follows from the particular form of d or d, see ©, that the 
surface =£? has genus zero and is thus homeomorphic to a sphere [17] . As a consequence this 
Riemann surface can be rationally uniformized, in the sense that it is possible to find two 
rational functions, say tt and tt, such that J2 = {(tt(s),tt(s)) : s G C U {oo}}. Furthermore, 
as shown in Chapter 6 of [13], we can take tt(s) = [x4 + xi]/2 + ([x4 — xi]/4)(s + l/s), x\ and 
X4 being defined below (jSJ); it is then possible to deduce a correct expression for tt, since by 
construction the equality Q(tt,tt) = has to hold. For more details about the construction 
of Riemann surfaces, see for instance [17j . 



2.2 Uniformization and meromorphic continuation 

But rather than the uniformization (tt, tt) proposed in [13] and recalled at the end of the 
previous subsection, we prefer using another, that will turn out to be quite more convenient. 
This new uniformization, that we call (x, y), is just equal to (tt o L, tt o L), where 

L (z) = — , zq = — exp (—m/n) . 

z- z 

We notice that zq is such that tt(zq) = tt(zq) = tt(zq) = vr(^o) = 1 and that its explicit 
expression above is due to (H3), for more details see Remark 1111 Then, starting from the 
formulations of (tt, tt) and of L, we easily show that the expression of the new uniformization 
can be 

(z + z )(z + z^) „. / ^_( z + z o) 2 



x (z) 



(9) 



(z - Zq) (z - z ) ' {z - zo 

Compared to (tt,tt), this uniformization (x,y) has the significant advantage of transfor- 
ming the important cycles (i.e. the branch cuts [xi,X4] and [2/1,2/4], the unit circles {\x\ = 1} 
and {\y\ = 1}) into very simple cycles, since the following equalities hold, see also Figured) 

x" 1 ([xi,x 4 ]) = lU{oo}, x~ l ({\x\ = l}) = zMU{oo} 



y 1 ([yi,y4l) 

% In 
' 



z IRu{cx3}, y-\{\y\ = 1}) 



Z !lU {OO}. 



(10) 



1x1=1 



n 



13-1=1 



Figure 4: The uniformization space C U {00}, with on the left some important elements of 
it and on the right the corresponding elements through the coordinates x and y 

To obtain (|10p . it is sufficient to use the explicit expressions of the branch points xi,X4, 
2/1)2/4) see below (jSJ), as well as the explicit formulation of the uniformization, see ([9]). 

Let us go back to £ and rj, the automorphisms of the algebraic curve £2 introduced in 
Subsection II. li Thanks to the uniformization ([9]), they define two automorphisms ^ and 77 
on C U {00}, which are determined by 

C 2 1 c c x2 2 n Pi -iy 2 +Pi,oV 1 ni N 
£ =1, xo£ = x, yo£=— , rj =1, yor) = y, xorj= ■ . (11) 

y pi,oy + pi,-i x 
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Using the well-known characterization of the automorphisms of the Riemann sphere C U 
{oo}, Q and (fTT)) . we obtain that £ and r\ have the following expressions: 

£(z) = l/z, r](z) = exp(-2m/n)/z. (12) 

The expression above of rj in terms of n is due to the assumption (H3) , see Remark [TT] 
for more details. Note also that leading to these particularly nice analytic expressions of £ 
and r] is another very pleasant property of the uniformization (x,y). 

As in Subsection 11.11 we call the group generated by £ and r\ 

W n = (£,r ] ) 

the group of the random walk. In the context of this article, W n is isomorphic to the dihe- 
dral group of order 2n, i.e. to the group of symmetries of a regular polygon with n sides, £ 
and rj playing the role of the two reflections. 

We are now going to state and prove Proposition [9] which actually is the main result 
of Subsection 12.21 and that deals with the continuation of the generating functions h and h 
defined in ([5]). For this, we need to describe the action of the elements of W n on some cones 
of the plane as well as to find some fundamental domains of the plane for the action of W n — 
we say that D is a fundamental domain of the plane for the action of W n if U W £w n w (D) = C 
and if in addition the latter union is disjoint. 

Let us take the following notation: for 9\ ^ 2 , let 

A(0i,0 2 ) = {texp(i6) : ^ t ^ oo,0i ^ ^ 2 } 

be the cone with vertex at and opening angles 0i, 2 . In particular, A(9, 0) = exp(z0)M_|_ U 
{oo}. Thanks to (|12|) . we obtain that the action of £ and rj on these cones is simply given 
by f(A(0i,0 2 )) = A(-0 2 , -0i) and r ? (A(0i,0 2 )) = A(-0 2 - 2vr/n, -0 X - 2vr/n). These facts 
are illustrated on the left of Figure 
Define now, for k £ {0, . . . , n}, 

. . fk-l k \ / fc + 1 k \ 

D+ = A[ 7r,-7T , D k =A[ TT, TT . 

\ n n J \ n n J 

Sometimes, we will write D instead of Dq = Dq and D n instead of D+ = D~. Clearly, 

n— 1 n— 1 

Do U D n U |J D+ U (J = C U {oo}. (13) 

k=l k=l 

The definitions of and , as well as (fT3l) , are illustrated on the right of Figure EJ 

It is immediate that for any k G {1, . . . , n}, we have = ^,{D^_^) and = ^(D^^). 
In particular, for any 2k G {1, . . . , n}, 

K = ((€ V) k ) (Do) , D- k = (( V o k ) (Do) . 

Likewise, for any 2k + 1 G {1, . . . , n}, 

D t k+1 = (£ o {rj o £) fc ) (A) , D- = {rj o (£ o r?)*) (Z) ) . 

With (jX3j) , these equalities prove that U w& w„w{Dq) = C U {oo}, in such a way that 
Dq is a fundamental domain for the action of W n on C — this is not quite exact, since each 
half-line A(/c7r/n, kir/n), k £ {0, . . . , 2n — 1} appears twice in the union U w <=w n w (Do)- 

We are now able to state and prove Proposition [9l after the weak recall on the lifting of 
functions that follows: any function / of the variable x (resp. y) defined on some domain 
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Figure 5: Important cones of the uniformization space 

D C C can be lifted on {z G C U {00} : x(z) G D} (resp. {z G C U {00} : y(» G I?}) by 
setting F(z) = f(x{z)) (resp. F(z) = f(y(z))). 

In particular, we can lift the generating functions h and h considered in ([5|) and we set 
H(z) = h(x(z)) and H(z) = h(y(z)); they are well defined on {z G C U {00} : \x(z)\ ^ 1} 
and {z G C U {00} : |y(#)| ^ 1} respectively. Consequently, on {z G C U {00} : \x(z)\ ^ 
1, |y(z)| ^ 1} (which equals A(— 7r/2,7r/2 — n/n), see Figured]), the functional equation ([6]) 
yields -ff(z) + H{z) — x lo y :io {z) = — in the sequel, we will often write x lo y JO (z) instead of 
x(z) io y(z) jo . 

Proposition 9. The functions H{z) = h{x(z)) and Ft(z) = h(y{z)) can be meromorphically 
continued from respectively A(— ir/2, n/2) and A(— ir/2 — n/n, tt/2 — ir/n) up to respectively 
C \ A(tt, n) and C \ A(-7r — Tr/n,ir — n/n). Moreover, these continuations satisfy 

H(z) = H(£(z)), H(z) = H(n(z)), VzGC, (14) 

and 

( if zi A(7r-vr/n,7r) (15) 

H(z) + H(z)-x io yi°(z) = I x io yi° (w(z)) if z G A(vr - n/n, tt) (16) 

V w£W„ 

where l(w) is the length of w, i.e. the smallest r for which we can write w = w\ o • • • o w r , 
with W\, . . . , w T equal to £ or r/. 

Remark 10. As a consequence of Proposition the generating functions h and h can 
be continued as meromorphic functions from the unit disc up to C \ [1,^4] and C \ [1,2/4] 
respectively. Indeed, the formulas h(x) = H{z) if x(z) = x and h(y) = H(z) if y{z) = y 
define h and h not ambiguously, thanks to fi!4\ )- Moreover, since we have x(A(n, n)) = [1, X4] 
and y(A(n — n/n,n — n/n)) = [1,2/4], see \10\) and Figure^ the previous formulas yield 
meromorphic continuations on C \ [1,374] an d C \ [1,2/4] respectively. 

Proof of Proposition To prove Proposition [9] we shall heavily use the decomposition (|13p . 
and precisely, we are going to define H and H piecewise, by defining them on each of the 
2n domains D that appear in the decomposition (fl~3j) to be equal to some functions Hd 
and Ho; it will then suffice to show that the functions FL and H so-defined satisfy the 
conclusions of Proposition [9] 

* In Do C {z G C U {00} : \x(z)\ ^ 1, \y(z)\ ^ 1}, we are going to use the most natural 
way to define Hd and Hd , i.e. their power series; so for z G -Do we set Hd (z) = h(x(z)) 
and H Do (z) = h(y(z)). 
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Then, for k G {1, . . . , n — 1}, we define H D + , H D + on D~£ and H D - , -ffp- on D k by 



Vz € L>+ = £ : (z) = (£ (*)) , (z) = —H D + (z) + xV° (*) , 

Vz G = 7] (D+_i) ■ H (z) = H D+ ( V (z)) , H (z) = -H (z) + s V° (*) . 



* At last, for z € D n , we set H Dn (z) = H n - (£(2)) and H Dn (z) = H n + (rj(z)) 

n— 1 n— 1 



Therefore we have, for each of the 2n domains D of the decomposition ()13p . defined two 
functions and -ffj> Then, as said at the beginning of the proof, we set H(z) = Hjj(z) 
and H(z) = Hd(z) for all z G D and for all domains D that appear in (|13p . 

With this construction, (I14p and (I15p are immediately obtained. In order to prove (]16p . 
we can use the fact that it is possible to express all the functions Hp, Hp in terms of Ho , 
Hd and x^yi" only; we give, e.g., the expression of H D + , for any 2k G {1, ... ,n}: 



k-l k-l 



H Dt (z) = -H Do ((r? o £) k (z)) + x i0 V i0 (£ ° (V o P (z)) ~ £ (fa ° , (17) 



2fe 

p=0 p=l 

as well as that of H D - , for any 2fe G {1, . . . , n}: 

k-l k-l 
H D - k (z) = -H Do ((£ o r,) k (z)) + E xV (r? o (e o r?)P(z)) - £ x V° ((£ ° ^)P(z)) . (18) 

p=0 p=l 

As a consequence, we get Equation (fTB]) for even values of n. Indeed, for this it is actually 
sufficient first to add Hf) n (z) and H£> n (z), in other words the equalities (fT7|) and (|TS|) above 
for = n/2, then to notice that if n is even, (£ o r/) n / 2 = (77 o £) n//2 , next to use that for 
zeD n , H Do (^o V ) n / 2 (z)) +H Do (^or]) n / 2 (z)) = i^^^ot])"/ 2 ^)), see (US), and finally 
to notice that W n equals 

{ 1, ^, . . . , (r^)™/ 2 " 1 , Cr?, ■ ■ ■ , m n/2 ~\ ^Cr 72 - 1 ,»/,..., ^(^) n/2 ^ 1 , (Cr?) n/2 } . 
Likewise, we could write the expressions of 

H-D- ' H D + , , H D + , , , 

2k 2k 2k + l 2k + l 2k + l 2k + l 

in terms of He> , Ho an d x lo y JO and we would verify that Equation (|16p is still true for 
odd n. Proposition [9] is proven. □ 

Remark 11. We can now explain precisely why the assumption (H3) dealing with the values 
of the transition probabilities is both natural and necessary for our study. 

If we suppose (HI) but no more (H3), then the uniformization ([PjJ is the same, with 
zq = — [2pi i _i] 1 / 2 + i[2pi q] 1 / 2 . The transformations (QTjJ) of the important cycles through the 
uniformization are also still valid and the automorphism £ is yet again equal to £(z) = 1/z. 
As for 7], it equals r](z) = Zq/z; in particular, the group of the random walk W = (£,rj) is 
finite if and only if there exists an integer p such that z 2p = 1. In this case, if n denotes the 
smallest of these positive integers p, the group W is then of order In. 

If a such n does not exist, then there is no hope to find a fundamental domain for the 
action of the group W , neither to obtain any equality like U6\) . 

If a such n exists, then by using the fact that z 2 ™ 1 = 1, in other words the fact that 
(— [2pi i] 1 ' 2 + i[2pi j o] 1 ^ 2 ) 2n = 1, we immediately obtain that pi^ = sin(g-/r/n) 2 /2, for some 
integer q having a greatest common divisor with n equal to 1. 

In this last case, we have zq = — exp(— iqir/n); it is then easily proven that the domain 
bounded by the cycles x^Ql, X4]) and y -1 ([l, 2/4]), namely A(arg(zo), it) = A(ir — qir/n,7r), 
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is a fundamental domain for the action of W if and only if q = 1. In particular, having an 
equality like [T6\) is possible if and only if q = 1, see the proof of Proposition^ 

But it turns out that having an equality like lilb]) is essential in what follows, particularly 
in Section^ where we have to know very precisely the behavior of H(z) + H(z) — x l °y J0 (z) 
near and oo. In the general case, the formulations of H and H — as solutions of boundary 
value problems — are so complex that we are not able to pursue the analysis. 

For all these reasons, we assume here that p\ o = sin(7r/n) 2 /2 for some integer n, in 
other words nothing else but (H3). 

3 Harmonic functions 

Section [3] aims at introducing and studying a certain harmonic function associated with the 
process, which will be of the highest importance in the forthcoming Section HJ 

It turns out that this harmonic function will be obtained from the expansion near of 
^2weW„(~ ^Y^x lo y^°(w(z)), quantity which is appeared in (fT6j) ; this is why we begin here 
by studying closely the behavior of the latter sum in the neighborhood of 0. 

Note first that thanks to the expression (]12p of the automorphisms ^ and ij, we have 

n-l 

(-l) l(w) xV° (w (z)) = [x io y io (eM-^kir/n)z) - x io y io (exp(-2ikir/n)/z)] . 

WGW n k=o 

(19) 

Let us now take the following notations for the expansion at of the function x^y 30 : 

x i °y i °( z )=Y / Kp(io,jo)z p , (20) 

and notice that with ([9]), we obtain that for z close to 0, 

x io y i °(l/z)=J2^(ioJo)z p . (21) 

In a general setting, if / is holomorphic in a neighborhood of with expansion f{z) = 
E p >o/p^> then Efc=o /(exp(-2z/c7r/n)z) = Y%Zq f(exp(2ikir/n)z) = T ip >o n fn P z np - 
This is why, by using (f20l) and (|2~T|) . we obtain that the sum (fT9j) is equal to 



Y *V° (v> {*)) = E n (*o.io) -K^(io,h)] z np - (22) 



We are now going to be interested in the term corresponding to p = 1 in the sum (|22p . 
and we set 

, /. \ K n (io,j Q )- K n (io,j Q ) . . 

fn{io,Jo)=n — Tn • (23) 

(-1) I 

We are now going to prove Proposition [5l but before we study some of its consequences. 

Corollary 12. The Doob f n -transform process of (X,Y) never hits the boundary. 

Remark 13. The function f n defined in A23\) is quite explicit. Indeed, by using the Cauchy 
product, K n ( and therefore also f n via A23\) ) can be written in terms of the coefficients of the 
expansions of x andy atO, and these coefficients are easily calculated, see Equation $25\). We 
don't know the factorized expression of f n for general values of n. However, we have given 
in Subsection \1.3l as examples, the factorized forms 0//3, f<± and f$. 
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Remark 14. The explicit expression and the harmonicity of the function f^ have already 
been obtained by Biane in J^]/. 

The quantity f&(io,jo) also appears as a multiplicative factor in the asymptotic tail dis- 
tribution of the hitting time of the boundary of for the process (X, Y) associated with 
n = 4 and starting from the initial state (iq, jo). Indeed, in JI| \2T$ , denoting by t = inf{k ^ 
: X(k) = or Y(k) = 0}, it is proven that ^(i ,j )[r > k] ~ C f 4 (i , j )/k 2 , where C > 0. 

In particular, we can specify Corollary {T^ in the case n = 4. Indeed, using the following 
equality for I < k (obtained from the strong Markov property of the process (X, Y) ): 

P (ioJo) [(X (I) ,Y(l)) = (ij) \r>k} = F {Hhjo) [(X (I) , Y (0) = (i, j)} F p j) k ~J} , 

the asymptotic of [20] yields that the Doob f ^-transform process is equal in distribution 
to the limit as k — > oo of the process conditioned on the event [r > k] . 

Remark 15. Let n ^ 3. Proposition^ gives that there exists at least one positive harmonic 
function for the process (X, Y). Corollary [3] entails that f n is in fact the unique — up to the 
positive multiplicative constants — positive harmonic function for (X, Y). 

Proof of Proposition The fact that f n takes real values is immediate from its definition. 



For the rest of the proof of (i) we are going to use the following straightforward fact: for 
any f(z) = 1 + Ylp^i fp,lZ p , note 1 + Yj P >\ fp,i zP tne expansion at of f(z) % ; then f p>i is 
a polynomial of degree equal or less than p in i, with dominant term equal to ff ^i p /p\. In 
particular, f pi is of degree exactly p if and only if fi^ ^ 0. 

In our case, it is immediate from Q that K\(l, 0) = — 4cos(7r/n) ^ and ki(0, 1) = 
— 4exp(z7r/n) ^ 0. This is why, for any non-negative integer p, K p (i,0) is a polynomial of 
degree p in i; likewise, k p (0, j) is a polynomial of degree p in j. In particular, K n (i,j) = 
Sp=o K p(0> f) K n-p(i, 0) is a polynomial in i,j of degree n, with dominant term equal to 

a /ci(o,i) p p «i(i,or-p .. n _ p 

^ pi (n-p)\ 

In this way, we obtain that f n is a polynomial in i,j of degree n, with dominant term equal 
to (after simplification) 



22n+l n ~ 1 



y 1 p=l 



^ CI sm{jm/n) cos{ir / n) n ~ p f i n ~ p . (24) 



Assertion (i) follows then immediately. 

To prove |(ii)| it is enough to show that n n is harmonic. To show that, start by using 
the obvious equality x l ~ 1 y^~ 1 (z)Q(x(z),y(z)) = 0, which reads x % yi{z) = piflX^y 3 (z) + 
Pi.oa^VO) +Pi,-ix i+1 y j ~ 1 (z) +pi- 1 x i ~ 1 y j+1 (z). Then, (|20D yields that 

^2 i K p( i ^)-P^,o K p( i + 1 J)-Pi,o K p( i - 1 J)-Pi ,-iMi+l, j-l)-fi -i K p( i - 1 J+ 1 )} zP 

is identically zero; this means that all the k p , p ^ are harmonic, hence in particular K n . 

In order to prove (iii) we need to know explicitly the expansions of x and y at 0. From 
([9]), we immediately obtain these expansions: 

x(z) = l + ^— - V(-l) p sin(p7i7n)z p , y(z) = 1 + 4 V (-lfpexp (ipn/n) z p . (25) 

tan(7r/n) 

K 1 ' p^i p^i 



We show now the first part of (iii) , namely the fact that f n (i, 0) = for all non-negative 
integer i. As it can be noticed from (|25p . the coefficients of x are real. For this reason, for 
all integers i and p, K p (i,0) is also real and thus f p (i,0) = 0; in particular, f n (i,Q) = 0. 
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As for the second part of (iii) , namely the fact that for all j ^ 0, f n (0,j) = 0, we prove 
that K n (0,i) is real — however, it isn't true that for all j and p, n p (0,j) is real. 

In order to obtain n p (0, j) — that is, the pth coefficient of the Taylor series of y{z) 3 — we 
add all the terms of the form k Pi (0, 1)k P2 (0, 1) x • • • x K Pj (0, 1) with p\ + • • • + pj = p, this 
is nothing else but the Cauchy's product of the j series y(z). In other words, using (125p . we 
add terms of the form p\ x • • • x Pj(— exp(t[pi + • • • + pj]ir jn). As a consequence, 
K p (0,j) can be written as ip p {j){— l) p exp(ipTr / n) , with (p p (j) > if j > 0. 

In the particular case p = n, we obtain n n (0, j) = —ip n (j)(—l) n ; K n (0,j) is therefore real 
and, immediately, / n (0, j) = 0. 

We prove now (iv) With (|25l) . it is clear that the sequence kq(1,0), . . . , K ra _i(l,0) is 
alternating, in the sense that for all p 6 {0, . . . , n — 1}, (— 1) p k p (1, 0) > 0. In particular, it 
follows from general results on power series that the sequence Kq(i, 0), . . . , 0) is still 

alternating, for any i > 0. 

In addition, by using the Cauchy's product of x(z) 1 and y^z) 3 , we obtain that n n (i,j) = 
K n (i,0) + K n (0,j) + ' exp(ip-K/n)(p p (j)K n - p (i,0). Then, by definition of f n (hj) 

and by using the fact that K n (i,0) and n n (0, j) are real, we get 

n-l 

f n (i,j) = 2n(-l) n ^2 (-l) P sin(p7r/n)v3p(j)K n _ p (i,0). 

But we have already proven that (p p (j) > if j > and that (— l) ra_p K n _ p (i, 0) > if i > 0; 
above, f n is thus written as the sum of n — 1 positive terms, and is, therefore, positive. □ 

Proof of Proposition^ Recall that h(pexp(i9)) = p n s'm(n6), see Subsection 11.31 Setting 
u = pcos(6) and v = psm(8) then yields: 

(n-l)/2 

/ l ( u , u )= ^ C^+i(_i)V-(2p+iV P +i_ 

p=0 



Next we easily check that up to a multiplicative constant, h(4>(io, jo)) equals the dominant 

■\n- 



term of f n (i ,j ), namely (2 2n+1 /(n - 1)!) YT~{ C p n sin(pvr/n) cos(^/n)™-^- p , see flM}- 



□ 



Proof of Proposition Since / n is a polynomial of degree exactly n, see (i) of Proposition [5j 
it is clearly enough to prove that for n ^ 5, / n (2, 2)// n (l, 1) 7^ 2™. For this we shall find 
both / n (l, 1) and / n (2, 2) in terms of n, it will then be manifest that / n (2, 2)// n (l, 1) 7^ 2 n . 

Let us first prove that / n (l,l) = 8n 3 / tan(7r/n). The Cauchy product of x(z) by y(z) 
and the use of (|25j) yields 



n-l 



«n(l,l) = 4(-l) n+1 n+ [16(-l) n /tan(vr/n)] ^pexp(zp7r/n) sin((n - p)vr/n) 



P 



Using then the exponential expression of sin((n — p)ir/n) as well as the identity Yl P =i P x 
[nx n (x — 1) — x(x n — l)]/[(x — l) 2 ] applied to x = exp^7r/n) entails: 

K n (l, 1) = 4(-l) n+1 n + [8(-l) n /(ztan(vr/n))][-n(n - l)/2 + nexp(-z7r/n)/(2isin(vr/n))]. 

Using (j23|) finally gives the announced value of / n (l, 1). 

Thanks to calculations of the same kind than for f n (l,l), we reach the conclusion that 
/ n (2,2) = (16/3)n 3 [n 2 + 2 - 6/ tan(vr/n) 2 ]/[sin(7r/n) 2 tan(vr/n)]. Then it becomes obvious 
that for n ^ 5, / n (2, 2)// n (l, 1) 7^ 2 ra , which completes the proof of Proposition [6j □ 



14 



4 Asymptotic of the Green functions 



Sketch of the proof of Theorem^ We shall begin by expressing Gij as a double integral, 
using for this Cauchy's formulas and ©, see Equation (|26p . Then we will make the change of 
variable given by the uniformization ([9]) and we will apply the residue theorem; in this way, 
we will write Gij as the sum Grjj 1 + Gi,j,2 of two single integrals w.r.t. the uniformization 
variable but on two contours a priori different, see (|27p and (|28p . Then we will show, using 
Cauchy's theorem and Proposition [9l that it is possible to move these contours of integration 
until having the same contours for both integrals Gij t \ and Gij^- Finally, by using f)16|) we 
will obtain (j29[) . which is the most important explicit formulation of the Gij, starting from 
which we will get their asymptotic. In (|29p . Gij will be written as an integral on the contour 
exp(2#)M + U {oo}, for some 6 € [tt — vr/n, n]. 

After having chosen an appropriate value of 9 € [tt — ir/n,ir], see (I3ip . we will see that 
this is quite normal to decompose the contour into three parts, namely a neighborhood of 
0, one of oo and an intermediate part. Indeed, the function x{z) l y{z) 3 that appears in the 
integrand of (|29p is, on the contour exp(z#)lR + U {oo}, close to 1 near 0, oo and strictly larger 
than 1 elsewhere. Next, we will study successively these contributions in three paragraphs, 
using for this essentially the Laplace's method, what will conclude the proof of Theorem [TJ 

Beginning of the proof of Theorem [J\. Equation ((6|) yields immediately that the gene- 
rating function G of the Green functions is holomorphic in {(x,y) £ C 2 : \x\ < 1, \y\ < 1}. 
As a consequence and using again Equation ([6]), the Cauchy's formulas allow us to write its 
coefficients Gij as the following double integrals: 

o, - ^ /L, - ^ fL, hlx) l%L-.t vk **. («> 



[2m] 2 J J\ x \ =1 sty' [2m] 2 J J\ x \ =1 &V j Q(x,v) 

\y\=l \y\ = l 

where the circles {|x| = 1} = {\y\ = 1} = {exp(z0) : 6 £ [0, 27r[} are orientated according to 
the increasing values of 8. 

With (f26|) . we can thus write Gij as the sum Gij = Gij t i + Gij^, where 



Gi ' j:1 [2m] 2 I\ x \=i x* ly^yiQix^)^' 



i f h{ y ) r dx d i r i r dx 

[2m] 2 J\ y \=i V 3 y W= i^Q(x,y) V [2m] 2 J\ y \ =1 y j ~ jo J\ x \=i x^Q (x, y) ' 



We are now going to make in Gjji the change of variable x = x(z). For this, we notice 
that if A(7r/2,7r/2) = {it : t € [0, oo]} is orientated according to increasing values of t, then 
the equality x(A(tt/2, it/2)) = —{\x\ = 1} holds — in the sense of the orientated contours — , 
see (fTUj) and Figure HI In this way and using in addition the identity q(x(z)) = H(z), we 
get 

1 f H(z) f dy ,. . , 

Gi j i = / — —r^ / . — . — -—^ — -x (z)dz. 

[2m] 2 JA(n/2,n/2) x{z) 1 J lyl=1 yJ Q (x{z), y) 

But Q(x(z),y) = if and only if y e {y(z),x(z) 2 /y(z)}, see (fTT|h Moreover, if z belongs to 
A(7r/2, vr/2) \{0, oo}, then \y(z)\ > 1, see FigureHl The residue theorem at infinity therefore 
entails that for such z, j^ =1 dy /[y J Q(x(z),y)] = —2m/[y(z) J d y Q(x(z),y(z))]. Finally, we 
have proven that 

q. = _L f EM *M dz (27) 

hhl 2mJ Hn/2tn/2)x ( z ) i y( z yd y Q(x(z),y(z)) 

A similar reasoning yields 
Gl2 = -—[ H{z)-x(z) io y(z) j0 y'(z) ^ 

hh2 2m 7 A (-7r/2-7r/n -7r/2-7r/n) X (zf y (z) j d x Q{x(z),y(z)) 
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We are now going to explain why it is possible to move the contours of integration of 
both integrals (|27|) and (|28|) up to A(9,9), for any 9 € [tt — tt/ti,tt] — see Figure [6] below. 

Start by considering Giji in (|27p . Thanks to Cauchy's theorem, it is sufficient to show 
that the integrand of Gij t i is holomorphic inside of A(-7r/2,7r), domain which is horizontally 
hatched on Figured! let us thus prove this fact. 

On one hand, with (UJ and Q, on the domain A(jr/2, tt) we get x'(z)/d y Q(x(z),y(z)) = 
— ^/(2[pl I oPl,-l] 1 ^ 2 -^), and the latter function has manifestly no pole. On the other hand, 
it is possible to deduce from the proof of Proposition [9] that the only poles of H are at zq 
and ~zq. In particular, using ([9]), we obtain that for i or j large enough, H(z)/[x(zYy(zy] 
has no pole in A(-7r/2, tt). Therefore, for i or j large enough, the integrand of Gj^i has no 
pole in A(7r/2,7r) and we can thus move the contour from A(-7r/2,7r/2) to A(9,9), for any 
6 S [tt/2,tt\. Note that it isn't possible to move the contour beyond A(7r,7r), since A(7r,7r) is 
a singular curve for H — indeed, remember that A(tt, tt) = x^]) and see Proposition [9j 

By similar considerations, we obtain that it is possible to move the initial contour of 
integration of Gij,2 up to A(9,8), for any 9 € [vr — vr/n, 3ir/2 — ir/n]. 

In particular, if we wish to have the same contour of integration for Giji and Gij t 2, we 
can choose A(9, 9), for any 9 € [tt — vr/n, tt] = [tt/2, tt] n [tt — tt/u, 3tt/2 — ir/n]. 




= A(7C/2,7C) 

Mill A(%-n/n,3%/2-%/n) 

ffffl A(TC-7C/re,7C) 

- - A(0,0) 



Figure 6: Change of the contours of integration in the integrals (|27p and (128 



Using then the equality x' \z)/d y Q(x(z),y(z)) = —y'(z)/d x Q(x(z),y(z)), that comes 
from differentiating Q(x(z),y(z)) = 0, as well as (f2?]h ([28]) and (fT6|) — we can use (fT6|) since 
9 £ [tt — 7r/n,7r] and thus A(9,9) C A(ir — 7r/n,7r) — , we obtain the following final explicit 
formulation for Gij, 9 being any angle in [tt — vr/n, tt]: 



G, 



1 



47r[pi j0 pi -i] 1/2 Ja{0, 



M w ) ™»o, JO 



x*V° (to («)) 



1 



x(z) i y(2;) 



rdZ. 



(29) 



On the contour A(8, 9) C A(7r — ir/n, tt), the modulus of the function x(z) l y(z) J is larger 
than 1, see Figure [U Moreover, it goes to 1 when and only when z goes to or to oo. 
This is why it seems normal to decompose the contour A(9, 9) into a part near 0, another 
near oo and the residual part, and to think that the parts near and oo will lead to the 
asymptotic of Gij, while the residual part will lead to a negligible contribution. But how 
to find the best contour in order to achieve this idea? In other words, how to find the value 
of 9 £ [tt — tt/ti,tt] for which the calculation of the asymptotic of (|29|) will be the easiest? 

For this, we are going to consider with details the function x(z) l y(zy , or equivalently 



Xj/i(z) = ln(x(z)) + (j/i)ln(y(z)). 
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Incidentally this is why, from now on, we suppose that j/i € [0, M], for some M < oo. 
Indeed, the function Xj/i is manifestly not adapted to the values j/i going to oo; for such j/i, 
we will consider, later, the function {%/ j)Xj/i( z ) = (i/j) ^ n ( x ( z )) + ^ n (y( z ))- Nevertheless, 
M can be so large as wished and in what follows, we assume that some M > is fixed. 

With ([9]), we easily obtain the explicit expansion of Xj/i a t 0: 

Xj/i (*) = E ^ 0"/*) * 2p+1 , (j/i) = h 2p+1 + zt 2p+1 + 2 (j/i) ^ 2p+1 ] . 

^ (30) 

Likewise, again with we get that for z near oo, Xj/i( z ) = X^=o u 2p+i(j/i)^-/z 2p ■ 

Consider now the steepest descent path associated with Xjjh m other words the function 
Zj/i(t) defined by Xj/i( z j/i(t)) = t- By inverting the latter equality, we easily obtain that 
the half-line (1 / v\(j / i))R + U {oo} is tangent at and at oo to this steepest descent path. 

Let us now set 

Pj/i = 1M (3/i) = V [2 (z + z5 + 2 (j/i) zq)] . (31) 

With this notation, we now answer the question asked above, that dealt with the fact of 
finding the value of 9 for which the asymptotic of the Green functions (|29p will be the most 
easily calculated: we choose 9 = axg(pj/ i ) — note that, from the definition of zq and (f3"Tj). we 
immediately obtain that axg(pju) € [ir — 7r/n,7r] — , and the decomposition of the contour 
M0.0) is 

A(arg(/j j/i ),arg(/j j/i )) = {pj/i/\pj/i\) [0, e] U {pj/i/\pj/i\) ]e, l/e[ U (pj/i/\Pj/i\) [V e >°°]- 

According to this decomposition and to (|29jl . we consider now Grjj as the sum of three terms 
and we are going to study successively the contribution of each of these terms. 

Contribution of the neighborhood of 0. In order to evaluate the asymptotic of the 
integral (j29l) on the contour (pj/i/\pj/i\)[Q, e], we are going to use the expansion at of the 
function 

1 £ (-1)^)^0^0(^(2)). 

This is why we begin here by studying the asymptotic of the following integral: 

jdz, (32) 



pj/i/\Pi/i\)[°A x(zfy(zy 



k being some non-negative integer. Using the equality l/[x(z) l y(zy] = exp(— ixju(z)) as 
well as the expansion (130(1 of Xjli a t an d then making the change of variable z = Pj/it, we 
obtain that (|32p is equal to 

p k+i r /IPj/ll ^ ex p(_^)ex P l-iJ2^ P+ i(j/i) { Pj/l t) 2p+1 ]dt. (33) 

But thanks to (f30|) . \v2p+i(j/i)\ ^ 4(M + 1) and therefore for all t € [0, e/\pj/i\], we have 
*£^i ^ P+ l(iA)(Pj/^) 2p+1 | < ^ 3 4(M + 1)/(1 - e 2 ). This is why 



exp -i E ^2 P +i (j/i) (Pj/it) P+ =1 + (ie 3 ) , 

V ^ I 

the O above being independent of j/i € [0, M] and of t G [0, e/ |] • The integral (|33j) can 
thus be calculated as 

p k +*[l + 0(ie 3 )) J t k eM-it)dt=( Pj/l /i) k+1 [l + 0(ie 3 )] J t k exp(-t)dt. 
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In the sequel, we choose e = 1/i 3 / 4 , so that ie/\pj/j\ — > oo and 0(ie 3 ) = 0(l/i 5//4 ). 

One could be surprised by this choice of e; in fact, in the upcoming paragraph "Conclu- 
sion", we will see that in order to obtain the asymptotic of the Green functions along the 
paths of states (i,j) € ^\ such that j/i — > tan(7) e]0,oo[, it would have been sufficient to 
have 0(ie 3 ) = o(l), but for the paths € such that j/i — )■ 0, it is necessary to have 
0(ie 3 ) = o(l/i), what affords the choice e = 1/i 3 / 4 . 

Finally, we obtain that for this choice of e, the integral (132H is equal to 



' dz=( Pj/i /i) k+1 k\[l + 0(l/i 5 /% (34) 



Pi,i/\pj/i\)[o,e] x{z) l y{zy 

where the O is independent of j/i £ [0, M\. 

We are presently ready to obtain the asymptotic of the integral (|29p on the contour 
(pj/i/\pj/i\)[0,e\. First, in accordance with (|22j) . we have that this integral equals 

1 f z np-l 

j-j— — j^gn Kp «,,*) - *), / (pj/i/|pj/(|)M ^^j^- 

Thus clearly, with (|34p . we obtain that all the terms corresponding in the sum above 
to p 2 will be negligible w.r.t. the one associated with p = 1. In addition, by using the 
definition (|23p of the harmonic function / n as well as (|34p for k = pn — 1 and p ^ 1, we get 
that the integral (i29j) on the contour (Pj/j/IPj/iDfO, e] is equal to 

-1)" (n - !)!/„ (i , jb) * (Pi/iA)" [1 + 0(1A 5/4 )] • (35) 



Contribution of the neighborhood of oo. The part of the contour close to oo, namely 
(pj/i/\pj/i\)[l/e,oo], is related to the part (Pj/i/\pj/i\)[0,e\ via the transformation z i-> l/~z. 
Moreover, it is clear from ([9]) that for / = x, f = y or / = ^ we w n {—^) l ^x to yi°{w), 



f(l/z) = f(z). 

Therefore, the change of variable z i— > 1/z immediately gives us that the contribution of the 
integral (|29p near oo is the complex conjugate of its contribution near 0. 



Contribution of the intermediate part. Let A e be the annular domain {z G C : e ^ 
\z\ ^ 1/e}. According to Figured! for all z € A(ir — ir /n,ir)PiA 6 we have \x(z)\ > l + r] Xte and 
\y(z)\ > 1 + rjy <e , where rj Xje > and 7] v>e > 0. In fact, since x'(0) ^ and y'(0) ^ 0, we can 
take r]x,t > rje and r] y>e > r\e for some r\ > independent of e small enough. 
Let us now consider 



L = sup 

zG.A(tt— 7r/n,7r) 



E 



x l °y 30 (z) 



and let us show that L is finite. For this, it is enough to prove that the function s, defined by 
s ( z ) = [^2 w eW (—^) l( - w ^ x^y 30 (w(z))] / [x l °y 30 (z)] , has no pole in A(tt — ir/n, w) — including 
oo. But by using ([9]) and (fl~2l) . we see that the only poles of the numerator of s are the 
zo exp(2ip7r/n) for p £ {0, ... ,n — 1}. Among these n points, only zq is in A(ir — ir/n,ir). 
But in s, we have taken care of dividing by x lo y' JO {z), so that s is in fact holomorphic near 
zq. Moreover, it is easily shown that s is holomorphic at oo. Finally, we have proven that 
s has no pole in A(tt — n/n, 7r), hence s is bounded in A(tt — 7r/n,ir); in other words, L is 
finite. 
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The modulus of the contribution of (|29|) on the intermediate part (pj/i/\Pj/i\)]e, l/e[c 
A(tt — ir/n, 7r) Pi A € can therefore be bounded from above by 

1 1 L (36) 



3-30 



4vr[pi i oPi,-i] 1 / 2 e 2 (1 + r ? e) i " i o(l + r/e) 

Note that the presence of the term 1/e 2 in (|36|) is due to the following: one 1/e appears 
as an upper bound of the length of the contour, while the other 1/e comes from an upper 
bound of the modulus of the term 1/z present in the integrand of (129|) . 

Then as before we take e = 1/z 3 / 4 , and then we use the following straightforward upper 
bound, valid for i large enough: 1/[1 + 7//i 3 / 4 ]* ^ exp(— [r//2]i 1 / 4 ). We finally obtain that 
for i large enough, (f36|) is equal to 0(i 3 ^ 2 exp(— [r]/2]i l / A )). 

Conclusion. We have seen that the contribution of the integral ([29]) in the neighborhood 
of is given by ([35]) . that the contribution of (f29|) in the neighborhood of oo is equal to the 
complex conjugate of (I35p and that the contribution of the residual part can be written as 
0(i 3 / 2 exp(-[r//2]i 1 / 4 ). Therefore, with ([29]) and ([35]), we obtain that 

Gij = 1 — 1TN (-If in - l)!/ n (i ,h) i [(pj/i/iT - [Wi^T] +0(l/i n+5/4 ). (37) 

47r LPl,0Pl,— lj ' 



Moreover, starting from pip , we easily derive 

, ,.,n , M „ 2i(-l) n+1 sin(narctan( T iA-tan(vr/n))) 

[Pj/iM -ifii/ih) = & r ( i m, xr . u - 2 W2 • 

^ [cos^vr/nj^ [i A + 2ij) + j z \ 

The latter equality, ([3?]) and Remark [2] conclude the proof of Theorem Q] in the case of 

7 e [o,tt/2[. 

* Note that having o(l/i n ) instead of 0(l/i n+5 / 4 ) would have been sufficient for 7 G 
]0, 7r/2[, since in this case, Remark [2] implies that {pj/i/i) n — (~p~jji/i) n ~ K^ji n with 
if 7 / 0. 

* On the other hand, if 7 = then (pj/i/i) n — (fijji/i) n ~ K Q j/i n+1 with K ^ and 
it is necessary to have something like o(l/i ri+1 ) in (|37|) . as it is actually the case with 

0(l/i"+ 5 /4). 

To prove Theorem Q] in the case 7 = 7r/2, we would consider rather than 

and we would then use exactly the same analysis; we omit the details. 
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